FASTC-ID(RS)T-0549-92 


J-AZD  I  UO  I 


FOREIGN  AEROSPACE  SCIENCE  AND 
TECHNOLOGY  CENTER 


ELIMINATION  OF  RHYTHM  FOR  PERIODICALLY  CORRELATED  RANDOM  PROCESSES 


by 

K.S.  Voychishln,  Ya.  P.  Dragan 


Approved  for  public  release; 
Distribudon  unlimited. 


FASTC-ID(RS)T-05A9-92 


HUMAN  TRANSLATION 


FASTC-ID(RS)T-0549-92  14  January  1993 

ELIMINATION  OF  RHYTHM  FOR  PERIODICALLY  CORRELATED 
RANDOM  PROCESSES 


By;  K.S.  Voychishln,  Ya.  P.  Dragan 
English  pages:  7 


Source:  Otbor  i  Peredacha  Informatsii,  Nr.  33,  1972; 

pp.  12-16 


Country  of  origin:  USSR 

Translated  by:  Charles  T.  Ostertag,  Jr. 

Requester:  FASTC/TATE/lLt  Douglas  E.  Cool 

Approved  for  public  release;  Distribution  unlimited. 


DTIC  QUALITY  INSPECTED  3 


{  Aoo»s8ion  for  j  \ 

STIS  GRA4I 
DTIC  TAB 
Unannounoed 
JustiTlcatl( 

a 

□ 

in 

i 

Bv 

Distribution/ 

Availability  Codes 

Diat 

Avail 

Spec 

and/or 

ial 

THIS  TRANSLATION  IS  A  RENDITION  OF  THE  ORIGINAL 
FOREIGN  TEXT  WITHOUT  ANY  ANALYTICAL  OR  EDITO¬ 
RIAL  COMMENT  STATEMENTS  OR  THEORIES  ADVO¬ 
CATED  OR  IMPLIED  ARE  THOSE  OF  THE  SOURCE  AND 
DO  NOT  NECESSARILY  REFLECT  THE  POSITION  OR 
OPINION  OF  THE  FOREIGN  AEROSPACE  SCIENCE  AND 
TECHNOLOGY  CENTER. 


PREPARED  BY: 

TRANSLATION  DIVISION 
FOREIGN  AEROSPACE  SCIENCE  AND 
TECHNOLOGY  CENTER 
WPAFB,  OHIO 


1^93 


FASTC-ID(rs)T-0549-32 


14  January 


U.  $.  BOARD  ON  GEOGRAPHIC  NAMES  TRANSLITERATION  SYSTEM 


Block 

Italic 

Transliteration  . 

Block 

Italic 

Transliteration 

A  a 

A 

4 

A,  a 

P  P 

9 

R,  r 

B  6 

B 

4 

B,  b 

C  e 

C 

t 

S,  s 

B  ■ 

B 

4 

V,  V 

T  T 

T 

44 

T.  t 

r  r 

r 

t 

0,  g 

y  y 

y 

f 

U,  u 

A  A 

a 

9 

D,  d 

♦  ♦ 

P,  f 

E  • 

B 

4 

Ye,  ye;  E,  e» 

X  X 

X 

« 

Kh;  kh 

IK  M 

M 

X 

Zh,  zh 

U  tt 

u 

y 

Ts,  ts 

3  9 

s 

f 

Z,  z 

M  ** 

V 

V 

Ch,  ch 

H  H 

H 

y 

I,  1 

lU  a 

Ul 

itf 

Sh,  sh 

A  A 

ft 

i 

y,  y 

m  Q 

Hi 

Mf 

Shch,  shch 

K  K 

K 

K 

K,  k 

T>  % 

% 

% 

It 

J1  a 

n 

4 

L,  1 

U  M 

u 

If 

y,  y 

M  M 

M 

M 

M,  m 

b  fc 

b 

A 

1 

H  M 

H 

H 

N,  n 

3  9 

9 

> 

E,  e 

0  o 

0 

4 

0,  o 

10  b 

K> 

• 

yu,  yu 

n  n 

n 

n 

P,  P 

R  M 

ft 

4 

ya,  ya 

*yc  Initially,  after  vowels,  and  after):,  t» ;  e  elsewhere. 
When  written  as  8  in  Russian,  transliterate  as  y8  or  C. 

RUSSIAN  AND  ENGLISH  TRIGONOMETRIC  FUNCTIONS 


Russian 

Engl Ish 

Russian 

English 

Russian 

English 

sin 

sin 

sh 

slnh 

arc  sh 

slnh  * 

cos 

cos 

ch 

cosh 

arc  ch 

cosh_' 

tg 

tan 

th 

tanh 

arc  th 

tanh  ^ 

ctg 

cot 

cth 

coth 

arc  cth 

coth”* 

sec 

sec 

sch 

sech 

arc  sch 

sech  * 

cosec 

CSC 

csch 

csch 

arc  csch 

csch  * 

Russian 

English 

rot 

curl 

Ig 

log 

GRAPHICS  DISCLAIMER 

All  figures,  graphics,  tables,  equations,  etc. 
merged  into  this  translation  were  extracted 
from  the  best  quality  copy  available. 


1 


ELIMINATION  OF  RHYTHM  FOR  PERIODICALLY 
CORRELATED  RANDOM  PROCESSES 


Voychishin,  K.S.,  Dragan,  Ya.P., 
L  '  vov 

Submitted  10  June  1971 


After  a  Fourier  expansion  of  periodic  functions  into  harmonic 
components  it  ceased  to  be  considered  as  a  method  of  solution  of 
problems  of  mathematical  physics  and  it  was  given  a  physical 
interpretation  in  the  theory  of  oscillations,  the  physical  systems 
began  to  be  treated  as  sets  of  oscillators.  Then  came  the  time  of  a 
certain  f etishization  of  harmonic  oscillations,  the  role  of  which  grew 
all  the  more  in  connection  with  the  successful  application  of  the 
method  of  harmonic  expansion  to  the  investigation  of  periodic  movements 
of  celestial  bodies  and  with  the  development  of  radio  engineering. 

These  circumstances  served  to  a  considerable  degree  as  the  reason 
that  for  the  development  of  methods  of  analyzing  more  complex 
phenomena,  such  as  meteorological,  economical,  etc.  for  example, 
initially  the  same  ideas  of  expansion  of  functions  into  harmonics  were 
used.  Thus  emerged  periodic  grammanalysis ,  the  essence  of  which  is 
from  a  complex  oscillation  to  identify  the  simple  harmonics,  to 
determine  their  frequencies,  and  on  the  basis  of  this  to  judge  the 
presence  of  specific  oscillators  in  the  system  being  investigated. 

However,  it  was  shown  that  an  harmonic  expansion  has  a  direct 
meaning  only  for  linear  invariant  systems,  for  which  the  harmonics  are 
their  eigen- f unct ions .  Moreover  the  development  of  the  theory  of 
expansions  led  to  the  theory  of  stationary  random  processes  (SP),  which 
showed  the  inconsistency  of  periodic  grammanalysis  as  a  method  of 
studying  many  complex  natural  phenomena,  since  in  a  general  case  the 
harmonics  here  are  the  result  not  of  the  presence  of  any  oscillators, 
but  of  the  stability  of  conditions  of  flow  of  the  process  (statistical 
equilibrium) . 

As  a  general  model  of  a  "periodic"  (rhythmic)  phenomenon  one 
should  consider  [1]  periodically  correlated  random  processes  (PKSP), 
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the  statistical  characteristics  of  which  are  periodic.  Here  rhythm 
should  be  considered  not  as  a  result  of  the  presence  of  harmonic 
oscillators,  but  as  a  result  of  the  correlativity  of  SP  harmonics,  and 
the  period  of  rhythmical i ty  treated  as  a  period  of  this  correlativity. 
From  this  point  of  view  the  harmonic  composition  of  the  process 
indicates  only  the  possible  periods  of  rhythm.  For  example,  if  an  SP 
is  polyharmonic,  then  the  period  of  the  harmonics  has  to  be  a  multiple 
of  the  period  of  correlativity,  which  is  confirmed  in  particular  by  an 
analysis  of  the  reactions  of  systems  with  periodically  changing 
parameters  to  a  stationary  SP.  From  here  it  is  evident  that  to  find  a 
frequency  spectrum  for  the  purpose  of  exposing  the  rhythm  makes  sense 
only  for  a  polyharmonic  SP. 

At  the  present  time  the  statistics  of  SP  are  supported  only  by 
concepts  of  stationary  processes,  therefore  its  recommendations 
frequently  are  inadequate  for  investigations  which  are  nonstationary, 
rhythmic  phenomena  in  particular.  However,  the  need  arises  for 
reconsidering  these  recommendations  from  the  point  of  view  of  the 
theory  of  PKSP  as  a  model  of  a  rhythmic  phenomenon.  One  of  such 
recommendations  is  the  preliminary  separation  of  the  high-frequency 
components  of  SP.  For  eliminating  the  regular  periodic  trend  (cycle) 
from  the  random  process  S(t)  it  is  recommended  [5]  either  to  select 
the  recordings  of  its  values  in  the  same  point  of  the  cycle  over  the 
period,  or  to  take  for  analysis  the  average  variables  for  the  same 
period.  In  this  case  the  investigated  process  represents  an  additive 
model  of  a  random  rhythmic  phenomenon  ^(t)=  where  /(t)  - 

determinate  periodic  function,  and  (t)  -  stationary  process. 

The  purpose  of  this  work  is  a  detailed  investigation  of  the 
questions  brought  up  above  about  the  model  of  rhythmic  phenomena  in  the 
form  of  PKSP. 

Since  in  the  general  model  in  the  form  of  PKSP  the  rhythm  is 
included  in  the  correlativity  of  the  harmonics,  then  the  elimination  of 
rhythm  in  this  model  is  equivalent  to  obtaining  a  stationary  SP  from 
PKSP.  In  particular  it  is  known  [4]  that  the  readings  of  the  values  of 
PKSP  after  a  period  of  its  cor relati  /  ity  form  a  stationary  random 
sequence,  i.e.,  as  a  result  of  the  read-out  the  rhythm  is  eliminated. 

Actually,  if  J  (tQ+kT)  -  the  sequence  of  readings  {k=l,  A',  and  to 
-  initial  phase),  then  in  view  of  the  periodicity  of  mathematical 
expectation  of  PKSP,  i.e.,  the  validity  of  the  correlation 
m|  (t+T)=mj  (t)  we  obtain 


+  AT)  » 


(1) 
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the  sequence 


Since  the  covariation  of  PKSP  is  also  periodic,  i.e., 
bj(t  +  T,  u)=bj(t,  u),  then  for  the  covariation  of 

^  we  will  have 


(/*  +  An 5 (^0 + ^  -  *1  ('o + An  (A - /) n  -  iu.  (A  -  /in.  ( 2 ) 


It  is  evident  from  formulas  (1)  and  (2)  that  in  the  case  of  a 
fixed  to  the  random  sequence  {  |  (tQ  +  kT)  is  a  sequence  which  is 
stationary  in  a  wide  sense,  a  sequence  whose  characteristics  evidently 
depend  on  the  initial  phase  tQ. 


Averaging  the  PKSP  values  for  the  period  and  relating  these 
results  to  a  point  with  the  phase  tQ+kTm  we  obtain  the  sequence 
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which  can  be  interpreted  as  the  readings  over  the 


period  of  the  moving  average  with  averaging  for  the  period  of 
correlativity;  i.e.,  to  obtain  'yjlj=V|(to  +  kT).  In  this 

I 

case  the  characteristics  of  the  second  order  of  such  a  sequence  will 
be  the  mathematical  expectation  m'^  =">  q  {to  =  kT)  and  the  covariation  [6] 


r, (/+«./)  «  R^iu)  (  3  ) 

u-T 


where  Bq{ii)  -  null  correlation  component  of  PKSP  [1]. 

On  the  basis  of  the  resulting  formulas,  taking  into  account  that 
in  {t)=mo,  where  mo  -  null  coefficient  of  expansion  of  the  function 
m ( • )  into  a  Fourier  series  for  mathematical  expectation,  we  find 

=  m,. 


and  for  covariation  - 

(»-/+i)r 
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In  comparing  the  expressions  (1)  with  (4)  and  (2)  with  (5),  we  arrive 
at  the  conclusion  that  for  PKSP  the  selection  of  the  readings  after  a 
period  is  not  equivalent  to  the  selection  of  the  average  for  the 
period.  For  eliminating  the  dependence  of  the  averages  for  the  period 
on  the  initial  phase  it  can  be  assumed  that  there  is  averaging  for  the 
period  of  characteristics  (1)  and  (2)  with  respect  to  phase  tQ.  Then 
for  these  averages  we  obtain 


K'”l  =  "”0  ^  ^  ^ 

and  respectively 

r 

where  operator  of  averaging. 

0 

From  a  comparison  of  formulas  (4)  and  (6),  and  also  (5)  and  (7), 
it  is  evident  that  this  method  is  not  equivalent  to  the  selection  of 
the  averages  for  the  period,  since  only  the  mathematical  expectations 
are  equal.  However,  all  the  rhythm  methods  are  eliminated. 

There  are  specific  characteristics  in  the  application  of  these 
recommendations  (readings  over  a  period  or  averaging  for  the  period)  in 
the  case  when  the  null  correlation  component  BqC*)  is  periodic,  since 
its  period  Tq  should  be  a  multiple  to  the  period  of  cor r e 1 a t i v i t y , 
i.e.,  To=nT  (n  -  integer  [3]).  Since  formula  (4)  can  be  written  in  the 
form 

rS,  =  Y  j  ^0  “  T  j  ^0 

0  0 


then,  according  to  the  Bolyai-Bohr  theorem  [2],  in  the  case  when  the 
undetermined  integral  from  Bq  is  limited  it  is  a  periodic  function. 
Therefore  the  function  r as  the  difference  of  periodical  functions  is 
also  periodic. 


It  follows  from  the  properties  of  the  Fourier  expansion  that  the 
period  of  the  indeterminate  integral  of  the  periodic  function  is  the 
same  as  the  period  of  the  initial  function.  Therefore  the  discrete 
spectra  of  functions  Bq  { u )  and  R (u)  will  be  concentrated  on  the  same 


geometric  progression  of  frequencies  with  the  difference 
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since  the  readings  across  the  period  T  will  have  the  result  that  the 
spectra  of  sequences  Bq  ( RT )  and  R  (kT)  will  be  concentrated  on  the  set 
of  frequencies  {kAo»  k  =  0 ,  n  - 1  }  ,  then  it  is  evident  from  here  that  in 
this  case  the  averages  by  phases  and  the  averages  for  a  period  will 
produce  the  same  set  of  frequencies,  but  their  energies  are  different. 

For  an  illustration  of  this  from  formula  (3)  we  find  that  if 


5,(u) 


(8) 


then 


m 

s 


sin  mA^T 


Then,  taking 


into  account,  we  obtain 

2n 


'1-2  2 


sin  m- 


2n 


Keeping  in  mind  that  e‘^*^  =1  in  the  case  of  a  whole  q,  we  see  that  in 

the  last  expression  a  "sticking"  of  the  coefficients  of  expansion  takes 
place,  and  then,  designating 


we  obtain 
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For  the  covariation  of 
{ 7 )  we  find 


the  averaged  sequence  on  the  basis  of  formula 


.«  • 


or,  having  designated 
the  form 


we  convert  this  expression  to 
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(10) 
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A  comparison  of  formulas  (9)  and  (10)  confirms  what  was  said  above. 

The  meaning  of  the  last  correlations  becomes  understandable  when  it  is 
taken  into  account  that  if  the  null  correlation  component  is  periodic 
with  a  period  To=nT,  i.e.,  it  can  be  presented  in  the  form 

1(0-  i 


then  the  PKSP  will  be  polyharmonic.  The  correlation  function  of  such  a 
process 


bU.u)^  2 

or  after  simple  conversions 


where 

Since  all  B j ( • )  have  the  same  period  Tq,  as  hoi’),  then  b{t,  •)  as 
the  sum  of  the  periodic  functions  possesses  the  same  period  and  can  be 
presented  in  the  form 


where  functions  Ci^(')  are  periodic  with  a  period  Tq.  Therefore  the 
covariation  of  readings  across  the  period 

•”1  ifl 

*  c,(g. 

c,(g-  i 
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Finally,  for  averaging  by  phase  Iq  we  obtain 


r 

where  C)  -  y- J C, (/o) <«o- 

From  here  it  is  evident  that  with  any  of  the  methods  mentioned 
above  for  working  polyharmonic  PKSP  we  will  obtain  a  stationary 
approximation,  the  spectrum  of  which  is  concentrated  on  the  set 
{kAO'  k=0,  n-1}.  Then  as  a  result  of  the  correlation  To=nT  in  the 
given  case  it  is  possible  to  estimate  the  period  of  cor relat i vi ty . 

If  the  PKSP  is  not  po 1 yha rmonic ,  then  these  methods,  as  already 
mentioned  above,  give  a  stationary  approximation  which  does  not  contain 
information  about  the  rhythm  of  the  initial  process. 
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